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We rewrite the Lagrangian for a dilute Bose gas in terms of auxiliary fields related to the normal 
and anomalous condensate densities. We derive the loop expansion of the effective action in the 
composite-field propagators. The lowest-order auxiliary field (LOAF) theory is a conserving mean- 
field approximation consistent with the Goldstone theorem without some of the difficulties plaguing 
approximations such as the Hartree and Popov approximations. LOAF predicts a second-order 
phase transition. We give a set of Feynman rules for improving results to any order in the loop 
expansion in terms of composite-field propagators. We compare results of the LOAF approximation 
with those derived using the Popov approximation. LOAF allows us to explore the critical regime 
for all values of the coupling constant and we determine various parameters in the unitarity limit. 
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I. INTRODUCTION 

In 1911, Kamerlingh Onnes found that liquid ^He, 
when cooled below 2.2 K began to expand rather than 
contract [1]. The transition, later named the lambda- 
transition was recognized in 1938 as the onset of 
superfluidity [2, 3]. The connection with Bose- Einstein 
condensation (BEC), first argued by F. London on the 
basis of the near identical values of the lambda transition 
temperature and the critical temperature for BEC 
of noninteracting bosons[4, 5] sparked a series of weakly 
interacting BEC studies when Bogoliubov[6] pointed out 
that the BEC elementary excitations satisfy the Landau 
criterion for superfiuidity[7]. In the weakly interacting 
limit, the interactions can be characterized by a single 
parameter [8] — the scattering length a — giving the re- 
sults a powerful, general applicability. The hope of study- 
ing bosons with short-range inter-particle interactions of 
a strength that can be tuned all the way from weakly 
interacting {p^/^a <C 1) to universality (p^/'^a ^ 1), ap- 
peared thwarted when it was found that the three-body 
loss-rate in cold atom traps scales as oc near a Fesh- 
bach resonance[9, 10]. In cold atom traps, only fermions 
have been obtained in the strongly interacting, quan- 
tum degenerate regime in equilibrium[ll], in which case 
three-body loss is reduced by virtue of the Pauli exclusion 
principle. Recently, however, it was pointed out[12] that 
three-body losses can be strongly suppressed in an opti- 
cal lattice when the average number of bosons per lattice 
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site is two or less. The development of novel cold atom 
technology[13, 14] leads to the prospect of studying finite 
temperature properties, such as the BEC transition tem- 
perature, Tc, superfluid to normal fluid ratio, depletion, 
and specific heat, at fixed particle density p. 

At finite temperature the description of BEC's even 
in the weakly interacting regime remains a challenge. 
Standard approximations such as the Hartree-Fock- 
Bogoliubov (HFB) and the Popov schemes, generally 
fall within the Hohenberg and Martin classification [15] 
of conserving and gapless approximations which imply 
that they either violate Goldstone's (or Hugenholz-Pines) 
theorem or general conservation laws[15]. Both these ap- 
proximations predict the BEC-transition to be a first- 
order transition, whereas we expect the transition to 
be second-order[16]. The calculation of Tc, first under- 
taken by Toyoda[17] to explain the difference between 
the lambda-transition temperature (2.2 K) and the 
r° (3.1 K) of the noninteracting BEC at the same den- 
sity, exemplifies the difficulties of understanding the the- 
ory near T^. whereas Toyoda found a Tc -decrease with 
increasing scattering length, K. Huang later pointed out 
that the calculation had a sign error, giving an increasing 
value of Tc [18]. However, Baym and collaborators[19, 20] 
noted that the Toyoda expansion involves an expansion 
in a large parameter. Their calculation found a linear in- 
crease of {Tc-T°)/T^ with pa?. The fact that the hefium 
lambda transition temperature falls below may be ex- 
plained by quantum Monte-Carlo calculations[19], which 
found that the critical temperature of a hard-sphere bo- 
son gas increases at low values of pa^, then turns over 
and drops below r° near pa? « 0.1. 

In this paper, we discuss in detail a theoretical descrip- 
tion that we introduced recently [21] to describe a large 
interval of p^/'^a values, satisfies Goldstone's theorem, 
yields a Bose-Einstein transition that is second-order. 
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gives the same critical temperature variation found in 
Refs. 19 and 20 but at a lower order of the calculation, 
while also predicting reasonable values for the depletion. 
This method then resolves many of the main challenges 
in describing boson physics over a large temperature and 
p^l^a regime and it's predictions will be available for ex- 
perimental testing in the near future. The approach we 
present here is different from other resummation schemes 
such as the large- A*" expansion (which is a special case of 
this expansion), in that it treats the normal and anoma- 
lous densities on an equal footing. 

In the following, we will discuss the general features 
that arise when rewriting the original theory in terms 
of composite fields. One aspect of this approach is that 
one can systematically calculate corrections to the mean- 
field results presented earlier[21] in a loop expansion in 
the composite-field propagators. We derive the Feynman 
rules for such an expansion using the propagators and 
vertices of the mean-field approximation. At each level 
of this loop expansion one maintains the features that 
the results are both gapless and conserving. The broken 
C/(l) symmetry Ward identities guarantee Goldstone's 
theorem order-by-order in the loop expansion in terms of 
auxiliary-field propagators [22]. 

In our auxiliary field formalism, we introduce two 
auxiliary fields related to the normal and anoma- 
lous densities by means of the Hubbard-Stratonovitch 
transformation [23, 24], utilizing methods discussed in the 
quantum field theory community [22, 25, 26]. This trans- 
formation has already been shown to be quite useful in 
discussing the properties of the BCS-BEC crossover in 
the analogous 4-fermi theory for the BCS phase [27-29]. 
The path integral formulation of the grand canonical par- 
tition function can be found in Negele and Orland [30]. 
The Hubbard Stratonovich transformation is used to re- 
place the original quartic interaction with an interaction 
quadratic in the original fields. An excellent review of 
previous use of path integral methods to study dilute 
Bose gases is found in the review article of Andersen[16]. 
The use of path integral methods to study various top- 
ics in dilute gases began with the work of Braaten and 
Nieto [31]. Path integral methods have recently been 
used to study static and dynamical properties of the di- 
lute Bose gases [29, 32-37]. An excellent summary of 
this approach and its connection to the more traditional 
Hamiltonian approach is to be found in the recent book 
by Calzetta and Hu [38] . We also point out that the 1 /N 
expansion, which is a special case of the method being 
proposed here, has a long history of use in high-energy 
and condensed matter physics [39, 40]. It has been used 
to calculate the critical temperature by Baym, Blaizot 
and Zinn- Justin [20]. This calculation gives the same 
result for as the method we are describing here. How- 
ever, our approach can be used at all temperatures. Cor- 
rections to the 1 /N result to calculating were obtained 
by Arnold and Tomasik [41]. 

The paper is organized as follows: In Sec. II we dis- 
cuss the auxiliary-field formalism and rewrite the La- 



grangian for weakly interacting Bosons in terms of two 
auxiliary fields. In Sec. HI we derive the loop expan- 
sion by performing the path integral over the original 
fields and then performing the resulting path integral 
over the auxiliary fields by stationary phase. In Sec. IV 
we find the leading-order loop expansion in the auxil- 
iary fields (LOAF) for the action. In Sec. V we set the 
auxiliary-field parameter d and discuss the leading-order 
effective potential for both the ground state and at finite 
temperature. In Sec. VI we discuss related mean-field 
approximations. In Sec. VII we discuss numerical results 
for the theory at finite temperature and varying dimcn- 
sionless coupling constant p^^^a. We compare the LOAF 
approximation to the Popov approximation in detail. We 
conclude in Sec. VIII. Finally, in App. A we discuss the 
connection between regularization of the effective poten- 
tial and renormalization of the parameters. In App. B 
we give the rules for determining all the Feynman graphs 
for the expansion using the mean-field propagators and 
vertices. 



II. THE AUXILIARY-FIELD FORMALISM 



The classical action S[ 4>, (/)* ] is given by 

5[</),r] = J[dx]£[c^,r], (1) 

where [dx] = dt d^x and where the Lagrangian density is 



ih 



<t>*{x) {dt^{x))-{dt4>*{x))<i>{x)] 



(2) 



<P*{x)[ 



2m 



fJ-o 



}0(x) 



Ao 



0(^)1' 



Here po and Aq are the bare (unrenormalized) chemical 
potential and contact interaction strength respectively. 
We introduce two auxiliary fields, a real field, x(3^)i 
and a complex field, A(x), by means of the Hubbard- 
Stratonovitch transformation [23, 24], utilizing methods 
discussed in Refs. 22, 25, and 26. In our case, the 
auxiliary-field Lagrangian density takes the form 



£^uM(l^\x,A,A*] = [x{x)-Xn coshe\(j){x)\ 



2An 



A{x) ~ Ao sinh6'(/)^(x) 



2l2 



(3) 



which we add to Eq. (2). Here 6^ is a parameter which pro- 
vides a mixing between the normal and anomalous den- 
sities. In Sec. VI, we will see that choosing = leads to 
the usual large- TV expansion which has only the auxiliary 
field X [25, 26]. In lowest order, = gives a gapless so- 
lution very similar to the free Bose gas in the condensed 
phase. If instead we choose 9 such that sinh0 = 1, then 
in the weak coupling limit our results agree with the Bo- 
goliubov theory [6, 16], which represents the leading-order 
low-density expansion. Of course all values of lead to 
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a complete resummation of the original theory in terms 
of different combinations of the composite fields. 
For an arbitrary parameter 0, the action is given by 



^[$, J] 



(4) 



+ / da; 



[da;] [da;'] cj)a{x)g-^\[x]{x,x') cj)\x') 
f Xi{x)x'{x) 



«(a;) J"(a;)}. 



with 



g-'\[x] = Six, x') { Go 1% + V\[x]{x) } , (5) 
h* - fioj 2m dt 



T/" r ^( \ —( x{x)cos^d — yl(a;)sinh( 

Here we have introduced a two-component notation using 
Roman indices a, 6, c, • • • for the fields 4>{x) and (jf (x) and 
currents j{x) and J*(a;), 

r{x)^{cj){x),<P*{x)) , <Pa{x)^{c^*{x),c^[x)) , (6a) 
3''{x) = {3{x),f{x)) , 3a{x)^{r{x),j{x)) , (6b) 

for a = 1,2, and a three-component notation using Ro- 
man indices i,j,k,--- for the fields x{x), ^(a;), and 
A*{x), 

X\x) = {x{x).A{x)/V2,A*{x)/V2) , (7) 
S\x) = {s{x),S{x)l^,S*{x)lV2) , 

and 

X^{x) = {x{x),-A*{x)/V%-A{x)/V2) , (8) 
S,{x) = {six),^S*{x)/V2,-S{x)/V2) , 

for i = 1,2,3. For convenience, we also define 
five-component fields with Greek indices $"(a;) = 
((/)"(a;),x*(a;)) and currents J" (a;) = ( j"(a:), S'*(x) ). 
These definitions define a metric rja^p for raising and low- 
ering indices. We use this notation throughout this pa- 
per. 

The action is invariant under a global U{1) trans- 
formation, 0(a;) e*"0(a;), A{x) — > e^'"A(x), and 
x{x) I— >■ x{x)- In components, the equations of motion 
are 

[ h — fiQ + x{x) coshf?] (f){x) — A{x) 0*(a;) sinhf? = j{x) , 
x(a;)/Ao = | (t>ix) [^ cosh6i - s{x) , 
A{x)/Xo = (/)2(a;) sinh 61 - S{x) . (9) 

We note that substituting x{x) and A{x) from the last 
two lines of Eqs. (9) (for zero currents) into the first line 
of Eqs. (9) gives the equation of motion for the field 0(a;) 
with no auxiliary fields [16]. 



Parametrizing the Green function Q as 

Gix, a;') = G*\x'^x')} ' ^^^^ 

and using 



[da;'] g-\x, x') g{x', x") = Six, x") , (11) 

we obtain the equations 

[/iQ -Ai + x(a;)coshe'] G(a;,a;') (12a) 
- A{x) K* {x,x') sinh e = 5{x, x') , 

[ho- fi + x{x)cosh0]K{x,x') (12b) 
-v4(a;)G*(a;,a;')sinh6' = 0, 

and the complex conjugates. Here, G(a;, x') and K{x, x') 
are the normal and anomalous correlation functions. 

III. AUXILIARY-FIELD LOOP EXPANSION 



The generating functional for connected graphs is 



Z[J] = e^'^M/'"' = A^y" D$ , 



,iS['S>,J]/h 



(13) 



with 5[<f>, J] given by Eq. (4). Average values of the fields 
are given by 



;$"(a;); 



h 1 



SZ[J] 



Z[J] SJaia 



6W[J] 



.1=0 



5Ja{x) 



.7=0 



(14) 



If we integrate out the auxiliary fields A{x) and x{x), we 
obtain the path integral for the original Lagrangian of 
Eq. (2). The strategy we will use here is to reverse the 
order of integration and first do the path integral over 
the fields 0°(a;) exactly and then perform the path inte- 
gration over the auxiliary fields by stationary phase to 
obtain a loop expansion in the auxiliary fields. Perform- 
ing the path integral over the fields 0°, we obtain 



Z[J] ^J^' j Dxe 



iS'otf[x,./]/(e?') 



(15) 



where the effective action is given by 
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^off[x, J] = 2 // [d^] [dx']U^)g[x]\{x,x')nx) 



[Ax][ 



Xi{x)x^{x) 
2Ao 
h 



+ X^{X)S\X) 



-Tr[ln[g-i(a;,a;)]]}. (16) 



Here x^{x) is defined in Eq. (7). As shown in Ref. 22, 
the dimensionless parameter e (which we eventually set 
equal to one) in Eq. (15) allows us to count loops for 
the auxiliary-field propagators in the effective theory in 



analogy with h which counts loops for the 0-propagator 
for the original Lagrangian. The stationary point Xo(^) 
of the effective action are defined by 6Seg[x, J]/Sxi{x) = 
0, i.e 



Ao 



{ \Mx)\^ +fiT^e{G{x,x)}/i} coshO - s{x) 
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{(j)l{x) + nK{x,x)/i} sinh6l-5(x) , (17) 



where 0o(a;) is given by 

ro[Xo,JKx)^ f[dx']g[xo]\{x,x')j\x'). (18) 



Both xo{x) a-nd Ao{x) include self consistent fluctuations 
and are functionals of all the currents J" (a;). Expanding 
the effective action about the stationary point, we find 



Sosix, J] = 5'off[xo, J 

x{x\x)-xh{^)){x'{x')^xi{^')) + 



(19) 



[dx] [dx']V^'[xo]{x,x') 



where V^,^ [xo] [x, x') is given by the second-order deriva- 



tives 



25.^'[Xo](x,x') = 



5x'-{x)5x^{x') 
1 

a; 



(20) 



Xo 



r)ij 5{x,x') + I{ij[xo\{x,x') 



evaluated at the stationary points. Here 11^ [xo](a;, x') is 
the polarization and is calculated in App. B. We perform 
the remaining gaussian path integral over the fields Xi by 
saddle point methods, obtaining 



eW[J]^So + S,i,[xo.J\ 
eh 



(21) 



2i 



dx]Tr[\n[V-\x^.J]{x.x)]\ 



where 5*0 is a normalization constant. From this we 
calculate the order e corrections to cj)"" — SW/Sja and 
X* = SW/SSi. Schematically, these one-point functions 
are shown in Fig. 1. The vertex function r[$] is con- 
structed by a Legendre tranformation (see for example 
Ref. 42) by 



r[$] = J [dx] Jaix) $"(a;) - W[J] 



(22) 



Here r[$] is the generator of the one-particle- irreducible 
(1-PI) graphs of the theory [43-45], with 



_5r[$]_ 



r{x) 



(23) 



(t)(0) (|)(0) 



(|,(0) (j,(0) 



— ^3^-(|)(0) + • — ^ 

? + 

"\(j)(0) 



(9(s2) 



/(B) = 



- (j)(0) 



(|,(0) 



+ (9(e2) 



FIG. 1. Feynman diagrams for (j> and x to first order in e. 
Solid and wavy lines correspond to the propagators of (j> and x. 
Dashed lines denote the zeroth-order (h'^'^K 



Keeping only the gaussian fiuctuations in M^[J], we find 
er[$] - \ jj[dx] [dx'] <Pa{x)Q'\xTb{x.x')<p\x') 

^j^^j { "^'^""l^y -^jr[Hg-\x]{x.x)]] 

;Tr[ln[2?-i[$]](a;,a:)]}+ero-f ••• , (24) 



2i 



which is the negative of the classical action plus self con- 
sistent one loop corrections in the 0° and x* propagators. 
Here, Fq is an adjustable constant used to set the min- 
imum of the effective potential to have finite reference 
energy. The effective potential VoS [$] is defined for static 
fields <i> by 



= - Vo + 2 0a V[x]\ -^^ - 

Tr[ln[g-i[x](a;,x)]] 
Tr[ln[p-i[<i>](x,x)]] + 0(e2), 



(25) 



2iVT 
eh 

2iVT 



where 



v[xTb 



XCOshO — ^ — Asinh^ 
—A* sinh 9 x cosh 9 — ^ 



(26) 



We will see below that for the static case, Q'^^[x\{x^x) 
and I?~^[$](a;,a;) are independent of x. 

For a system in equilibrium at temperature T, we Wick 
rotate the time variable to Euclidian time r according to 
the Matsubara prescription, t — > —ihr. Then the effec- 
tive potential becomes the grand potential Q[^] per unit 
volume, Vcff[$] rt[^]/V. (Details of the Matsubara 
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formalism can be found for example in Ref. 30.) So to 
leading order in e, the thermal effective potential is given 
by 



A ZAo 



(27) 



1 

2fW 



Tr[ln[g-i[x](2;,x)]], 



and where Vq is a normalization constant. At the next 
order we have the additional term 



Tr[ln[X'-i[$](a;,a;); 



(28) 



Here and throughout this section, we suppress the de- 
pendence of quantities on Q and the thermodynamic vari- 
ables ( T, /i, y ) . The thermodynamic effective potential 
Veff['i?o] is obtained by evaluating the effective potential 
at zero currents. From (23), this is when the fields $0 
satisfy 



<5Vcff[^o] 



= , for a = 1, • • • , 5. 



(29) 



We call these the "gap equations" in analogy with the 
corresponding equations in BCS theory. 

The Green functions are periodic with Matsubara fre- 
quency w„ — l-Knjfi with /3 — l/{kBT), and are ex- 
panded in a Fourier series. 



{r—r')—LUn i'^—'T'') ] 



(30) 



k,n 



Writing the Green function equation in k-n space as 

^-i[x](k,n)^[x](k,n) = l, (31) 

we find 



(32) 



Cfc + X cosh — iaj„ —AsmhO 

—A* sinh 9 + X cosh + iuj^ 



where ^fe = efc - A^o- So 

det[g-'[x]{Kn)]^ujl+u;l 

where 



= [6 + XCOsh6']2 - lApsinh^e , 



(33) 



(34) 



Stable solutions are possible for aj| > 0. The trace-log 
term then becomes 



^Tr[ln[g-Mx](.,x)]]=^5:in[ 



2 I 2 ■ 



k ^ 

f d^k ( LUk 1 



} 



{f^Ml-e-l} 



(27r)3 I 2 P 



(35) 



So from (27) the effective potential to leading order in 
the auxiliary field loop expansion (LOAF) is given by 



Vefr[$] = Vo + [xcosh0-fio[ 
- ^ [(1)*^ A + (t)^ A*] sinh 6* 



(36) 



x' - 1^1 



<^ — + - In 1 - p"^"" 

2Ao J (2^)3 1 2 /3 ^ 



It is useful to introduce new variables x' smd A' as 

x' = X cosh 6 — fiQ , and A! = A sinh . (37) 
Then the effective potential (36) becomes 

1 



Veff[$'] = Vo + |<^Px'■ 
(x'+Mo)' , 



2 



2Ao cosh^ 6* 2Ao sinh^ d 



(38) 



(27r)3 



where now a;| = ( ffc + x' )^ ~ l^'P- The gap equations 
(29) are now written as 



Xo + Mo 
Ao cosh'^ 6 
A' 



2 , / d^fc efc-^xo 



(27r)3 2wfc 



= , (39) 

27l(/?Wfc) + l] , 



An sinh Q 



A' 



d3/c [2n(/3a;fe) + 1] 



(2vr) 



2wfc 



where 7i(a;) = l/(e^ — 1) is the Bose-Einstein particle 
distribution. The solutions <i>Q" of Eqs. (39) substituted 
into Eq. (38) determine the effective potential. 

To calculate the finite temperature effective poten- 
tial to order e we need to determine Tr[ln[2?~^(a;, a;') ] ]. 
For the static case in the imaginary time formalism, 
'Dij\x\(x,x') and nij[$](a;, x') are expanded in Fourier 
series' analogous to Eq. (30). So from Eq. (20) we obtain 



P,j[$](k,n) = ^+n,^-[$](k,n) 

-^0 



(40) 



IV. THE EFFECTIVE POTENTIAL IN THE 
CONDENSATE PHASE TO LEADING ORDER 

In the language of broken symmetry, the condensate 
phase is a phase where the U(l) symmetry of the theory is 
broken since then ( </>) ^0. From Eq. (38), the minimum 
of the effective potential is when 



X 00 



A', 



0. 



(41) 



Because of the C/(l) gauge symmetry, we can choose (/)o 
to be real, which means that A is also real. Hence, we 
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have the broken symmetry condition x' = A', and the 
dispersion relation reads 

Lol - Efe (e, + 2A') , (42) 

The latter is a consequence of the the Hugenholz-Pines 
theorem which assures that the dispersion relation does 
not exhibit a gap. This is equivalent to the Goldstone 
theorem for a dilute Bose gas with a spontaneously- 
broken continues symmetry. This connection is discussed 
in detail in Ref. 16. In the absence of quantum fluctua- 
tions in x' = A', one obtains the Bogoliubov dispersion, 
ujI — efc(e/c + 2A(/)q), by setting A' — Xcj)^ sinh^ ^ and 
sinh 6 = 1. 

In the spontaneously broken phase, the effective po- 
tential is 



Vefftx'] = Vo 



X 



12 



2An cosh^ 9 2An sinh^ 6 



(43) 



-t- 



2 



1 



Infl 



(27r)3 L 2 ' /3 
where x' is determined by the equation 



dx' 



d^fc 



Ag sinh 9 Ao cosh 9 



(44) 



(27r)3 



[2n{l3ujk) + I] =0 



These equations for Veff[x'] and x' contain infinite terms 
that need to be regulated. In order to regulate the ef- 
fective potential, we first expand ujk in a Laurent series 
in Cfe 



= V(efe + x')'-|^T = efe + x'-^ + 



'|2 



2efc 



(45) 



around fc — >■ (X). The first three terms in the series are 
responsible for the divergences in the integral in Eq. (38). 
To regularize the theory, we subtracting these three terms 
from uj]^ in the integrand, and replace the constant Vq 
the bare interaction strength Aq and chemical potential 
/j,o by regulated ones. This procedure gives the regulated 
effective potential 



V«f[$']=VR + |0|2x' 
(X' + A*)' , 



\Ar 



,* 2 



A' 



2A cosh^ 9 
d^fc r 1 
(27r)3 I 2 



2A sinh^ 9 
\A'? 



w/c - X + 



2efe 



/3 



ln[l 



(46) 



which is now finite. Similarly, the regulated gap equa- 
tions for A' and x' £^re now give as 



X + M 
A cosh^ 9 

A' 
A sinh^ 9 



r d^k r 
7 (2^1 



2wfe 

d^fc (2n{l3uJk) 



J (27r)3 I 



2n(/?a;fe) + 1 
1 



2wfc 



-} (47) 



which are also finite. 

This regularization scheme is equivalent to dimensional 
regularization as done for example in Ref. 46, or to con- 
ventional renormalization of the coupling constant and 
chemical potential as described in the review article of 
Andersen and discussed in detail for the LOAF approxi- 
mation in the App. A. 



V. SETTING THE PARAMETER 9 

Up to this point, we considered a one-parameter class 
of mean-field approximations governed by the parame- 
ter 9. The dispersion relation in the condensate phase 
for the leading-order auxiliary field (LOAF) approxima- 
tion is given by Eq. (42) 



; = v/e/c(efc + 2A sinhf 



(48) 



where A/X = {(j)'^) smh9 = [0^ + hK{x,x)/i] smh9. 
Now, we will choose 9 by demanding that in the weak 
coupling limit, when K{x,x) can be ignored, the disper- 
sion relation agrees with the one-loop low-density result 
obtained by Bogoliubov. Using a Hamiltonian formalism, 
Bogoliubov assumed 



4> = (pQ+ij] . 



(49) 



subject to the constraint ( -0 ) =0. Realizing that 00 ~ 
-\/]V, he then wrote the theory in terms of the classical 
Hamiltonian plus a quadratic fluctuation Hamiltonian, 
which he diagonalized. Using Eq. (49) and limiting to at 
most quadratic fluctuations, one has 

+ v*)(0o + ^)]' ^ {KM^ 4 ^S-^o ) 

+ V'V(0S</'o) + 0S^o(0S0o). (50) 

The minimum of the classical Hamiltonian defines /i = 
A (0Q0o)- One can reformulate[16] the Bogohubov theory 
in path integral language as the classical approximation 
plus gaussian fiuctuations. The inverse Green function 
in the gaussian fluctuation approximation now has 

where V""b[0] is defined in Eq. (5). This leads to the 
dispersion relation at the minimum: 



Wfe = \Jf-k{tk + 2A0^) 



(52) 



We will choose 9 such that our result for ojk reduces to 
the Bogoliubov dispersion relation (52) when we ignore 
quantum fluctuations in the anomalous density. This sets 
sinh 6* = 1 and coshfl = \/2. With our choice of 9, the 
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renormalized effective potential can be written as 



4A 



2A 



d^fc f 1 
(27r)3 I 



uJk- ek-X + — 



ln[l-e" 



(53) 



where now x' = X ~ M ^-^^d 



^^ = (£.+x')'-|^r (54) 

The equations for the auxihary fields are obtained from 
<5Vfff[$]/<5xJ-0, as 



A 
A 

2A 



^2 , d3fc ^ [2n(/3c^fc) + 1] 1 
^'(2^1 ^ 



2wfc 



2eJ 



(55b) 



From Eq. (41) we know that at the minimum of the ef- 
fective potential we have (x' — A) (j)o = 0, and we can 
replace /i by the physical density using 



P = 



2A 



(56) 



In the broken symmetry phase we have x' = ^ in which 
case Eqs. (55) become 



X' _^ , /'d^fc f [2n(/3^fc) + 1] 1 T 

^-^° + /(£F{^[2"(^'^^) + i]4}'(^^b) 

where po = 0o ^^^^ condensate density. 



VI. RELATED MEAN-FIELD 
APPROXIMATIONS 

For comparison, we will review next two related mean- 
field approximations. We will focus on the leading-order 
large- TV approximations, which corresponds to the choice 
of = in our formalism, and the Popov approximation 
that is widely used in the study of BEC condensates. 



A. Large-A'^ approximation in leading order 

The large-iV approximation corresponds to the value 
9 = 0. In obtaining the large-N approximation, one 
rewrites 4>*4> in terms of two real components and ex- 
tends the theory to N real components. The 0{2) [U{1)] 



symmetry is then extended to 0{N). Here the compos- 
ite field is X = ^4'i4>i/^- With appropriate rescaling, 
one can show[22] that the composite-field propagator is 
proportional to 1/N, so counting loops of bound-state 
propagators yields the 1/N expansion. In lowest order 
we find that this approximation in the BEC phase leads 
to the free-field dispersion relation. A related large- iV ex- 
pansion for the Bose gas at the critical temperature has 
been used successfully to characterize the behavior near 
the critical point [20]. One simplicity of this expansion is 
that the noninteracting-like dispersion relation simplifies 
the integrals present in the theory, and one can obtain 
analytic results even at finite temperatures. As with the 
general 9 result, the large- TV expansion also provides a 
complete resummation of the original theory. 

The large-iV finite-temperature effective potential in 
leading order is given by 



(55a) Vln[$]=Vo + xI 



^-^^-iTriHG-^]]. (58) 



The Matsubara inverse propagator in momentum space 
is now diagonal: 



lUJn - Efe - X 







-iuJn - Cfe - X 



We write the temperature-dependent last term in 
Eq. (58) as 



1 



Tr[ln[g- 



WJ (2^ ,^ 



(2^) 



/3 



ln[l - 



(59) 



where ujk = e/c + X- Inserting this into Eq. (58), the 
effective potential for the large- case is given by 



Vln[*] = Vo+xl 



,2 ix + p-oY" 



2Ao 

^ d^fc f LUk 1 

(27r)3 I T + /? 



(60) 



ln[l - 



Setting the derivative of the effective potential with re- 
spect to X equal to zero yields the gap equation. 



Ao 



A ^3 



(Pk 2nil3ujk) + l 



(2^) 



(61) 



The large-A^ effective potential in leading order 
is renormalized following the procedure discussed in 
Ref. [16]. We recognize that the infinite constant is re- 
lated to the renormalization of the chemical potential, 
i.e 



Ao " A 



A n3 



d^fc 1 



(27r)3 2 



(62) 



The renormalization is a consequence of the lack of a 
normal-ordering step in the the path-integral formalism. 
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in contrast with the usual Hamihonian formahsm. Per- 
forming the renormahzation, one obtains the finite gap 
equation, 



X + M 
A 



(2^ 



n{l3uJk) 



(63) 



where we have set the condensate density po = 10 p. 
Eq. (63) determines %[(/)] imphcitly, which is then re- 
inserted into the expression of Vln, so that Vln[^] be- 
comes solely a function of po — \(j>\'^- The renormalized 
potential is now 



(X + A^)' 



2A 

d^fc 1 



ln[ 1 - 



(27r)3 13 

The minimum of the effective potential is when 
9 Vln 



0, 



4>x = 0- 



(64) 



(65) 



So, in the large- A'^ mean-field approximation, the broken- 
symmetry regime, 7^ 0, corresponds to the condition 



(66) 



which gives the dispersion relation, uj^ — ^k, which is the 
same as the free-field theory dispersion. 

At finite temperature, the gap equation at the mini- 
mum 



^ Po = M - A 



(27r)3 

which gives the chemical potential as 



n{(3uJk) 



M = A { PO + ^ T'/' C(3/2) } 



Correspondingly, the phase transition 
at the free-field critical temperature: 



2fi 



AC(3/2)y^ 



2/3 



(67) 
(68) 

0) takes place 
(69) 



At the minimum, % = so that the value of the effective 
potential at the minimum as a function of temperature 
for T < is 



V, 



LN 



"2A 



d^fc 1 



(27r)3/3l^^[l^ 
ll r3/2 V^C(5/2) 
2A 



(70) 



167r2 2 

The total number density is determined from 



9 Vln 

dp. 



A 



(71) 



Hence, by combining Eqs. (68) and (71), we obtain the 
density of particles in the condensate as 



Po = P 



- r3/2 c(3/2) . 



(72) 



In summary, below Tc the large- iV approximation gives 
essentially the same results as a non-interacting gas since 
X = 0. Above Tc, the large-A'^ approximation gives rise 
to a self-consistent correction to the dispersion relation. 
The large- iV result above Tc is the same as that of the 
Popov approximation we review below. We will also find 
that the large- TV approximation is equal to the LOAF ap- 
proximation we are proposing here at high temperatures 
in the regime where ^ = 0. 



B. Hartree and Popov approximations 

The Hartree approximation is a truncation scheme that 
ignores correlation functions beyond the first two. Tech- 
nically this is obtained by setting the third derivative of 
the generating functional of connected graphs with re- 
spect to the external currents to zero, i.e. 



Sj{x)Sj{y)5j{z) 



= 



(73) 



Then, the vacuum expectation value of the expectation 
value of the field (/)[j](x) in the presence of external 
sources is 



{h-fi)(l>{x) + Xo\cf>(x)\^cl>{x) 
+ 2XohG{x,x) cj){x)/i 
+ XQhK{x, x) (j)*{x)li — j{x), 



(74) 



where h was defined in Eq. (5). Here <j){x)^ G{x,x') and 
K{x,x') are considered functionals of the current j{x). 
G{x, x') and K{x, x') have the same meaning as the nor- 
mal and anomalous correlation functions in Eq. (10). In- 
troducing new auxiliary fields xi^) ^-^d ^i^) by the def- 
inition. 



Xjx) + Pq 
2Ao 

Ao 



:^\(l>{x)\^ + hG(x,x)/i, 
= [(j){x)]^ + hK{x,x)/i . 



(75a) 
(75b) 



and setting j{x) = in Eq. (74) gives an equation for the 
average field, 

[ h + x{x) - 2A I 4>ix) |2 ] (t>ix) + A{x) r {x)^0. 

and its complex conjugate. Functional differentiation of 
Eq. (74) with respect to j{x') and j*{x'), ignoring third- 
order functional derivatives leads to equations for the 
Green functions G{x,x') and K{x,x'). We find 

[h + x{x)] G{x,x')+A{x)K*{x,x') ^ d{x,x') , (76a) 
[h + x{x)] K{x, x') -f A{x) G* (a;, x') = , (76b) 
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and the complex conjugates. Eqs. (76) can be written in 
matrix form as 



dx Q ^{x,x')Q{x' ,x") = S{x,x") 



(77) 



where 



g-(x,x')=^(x,x')(' + fif . (78a) 

The renormahzed eiTective potential for the Hartree ap- 
proximation can be written as: 



4A 



(79) 



1 

2A" ^ 2 
d^k r 1 r 
(2^ I 2 



'A* 



1 



/3 



Infl 



Wfc - Efc - X 

The minimum of the potential is given by 

Xo4>o-2\\M^(f>o + Ao(f>*o=0. (80) 



}■ 



d(f>* 



00 



Again, the U{1) gauge symmetry, allows us to choose (f)Q 
to be real at the minimum. Then according to Eq. (80), 
Aq is also real. Hence, Eq. (80) becomes 

[xo-2A|0oP + ^o]'/'o = O. (81) 

In the broken symmetry case, 4>o 7^ 0, we have 

Xo+Ao = 2X\cPo\\ (82) 

and the dispersion relation is 

u^l + 2A l^oP ] [e, - 2 ( Ao - A |0oP ) ] ■ (83) 

The Hartree approximation has the defect of not being 
gapless. This can be fixed by hand by ignoring the fluc- 
tuation in the anomalous density, that is by arbitrarily 
setting 



Ao-X\<f>o\' 



0. 



(84) 



This further approximation is known as the "gapless" 
Popov approximation [47]. 

The Popov approximation includes the self-consistent 
fluctuations of x, but treats A classically. Below Tc, the 
Popov approximation has the dispersion relation 



= <^-k{(k + 2Apo ) , 



(85) 



and the chemical potential is = 2Ap — Apo- The con- 
densate density pq is given by 



In the Popov approximation the critical temperature is 
the same as in the free- field case, Tc = Tq. 

The Popov approximation is the most commonly used 
mean-field theory of weakly interacting bosons at fi- 
nite temperatures[16]. This approximation has a gapless 
spectrum but is known to produce an artificial first-order 
phase transition as we shall see below. Unlike the LOAF 
and the Hartree approximations, the equations for 0, A, x 
in the Popov approximation are not derivable from an ef- 
fective action. 



VII. MEAN-FIELD RESULTS AND 
DISCUSSIONS 

We begin by comparing the predictions of the LOAF 
and Popov approximations for zero-temperature condi- 
tions. In the broken-symmetry phase at T = 0, we have 
x' = A, and the effective potential is given by 

1 + 



Vcff[x'] 



2A 



X 



(87) 



1 

47r2 



dk 



12 



1 

2A 



(x' + m)' ,2 



2ek 

2^2 5/2 
157r2 ^ 



Setting dVctf/dx' = 0, then gives x' as a function of fi at 
the minimum. We have 

X' = - A 1^ X"/' ■ (88) 

The above cubic equation can be solved explicitly. In 
particular, in the weak coupling limit we obtain 



x' = M - A 



2^2 
3^' 



3/2 



Using Eq. (56), in weak coupling we obtain 



P = 



A 



1 



37r2 



(89) 



(90) 



which agrees with the one-loop result corresponding to 
the original Bogoliubov approximation (see e.g. Eq. 89 
in Ref. [16]). By inverting Eq. (90), we derive fi{p) at 
weak coupling, as 



1 



1 \p 
3^VT 
32 



Ap 

= 8npa(^l + y ^/pa?J^T^ , 



(91) 



where we have set A = Anh? a/ra with a the s-wave scat- 
tering length. 

We can also calculate the condensate depletion, defined 
as p — Pq- From Eq. (57b) at T = 0, we obtain the exact 
LOAF result 



1 



P- Po 



X 



'3/2 



(92) 



10 




10^ 10^ 



FIG. 2. (Color online) Comparison of the predictions of 
the Popov and LOAF approximations regarding the zero- 
temperature values of the densities, x' = A, condensate frac- 
tion, po/p, chemical potential, p., and effective potential, Vcb, 
as a function of dimensionless parameter, p^^^a. In the case 
of the LOAF approximation, the normal and anomalous den- 
sities are equal, x' ~ whereas in the Popov approximation 
we have A = Apo. Note that there is no effective potential in 
the Popov approximation, because this approximation is not 
derivable from an action. 



Hence, using Eqs. (89) and (91) we obtain the weak- 
coupling result for the fractional depletion (see e.g. Eq. 
22.14 in Ref. 48) 




(93) 



first obtained by Bogoliubov in 1947 [6]. 

In Fig. 2 we depict the coupling constant dependence 
of the zero-temperature values of the normal densities, 
x', condensate fraction, pa/p, chemical potential, /i, and 
effective potential, Vcs^q. The coupling constant depends 
linearly of the dimensionless parameter, p^/^a. We note 
that in the case of the LOAF approximation, the normal 
and anomalous densities are equal, x' = whereas in 
the Popov approximation we have A = XpQ. Also, in the 
Popov approximation there is no effective potential, be- 
cause this approximation is not derivable from an action. 




FIG. 3. (Color online) Normal density, x') and anomalous 
density. A, from the LOAF and Popov approximations. The 
comparison between the LOAF and Popov approximations is 
carried out for p^^^a = 1, p^^^a = 0.4, and p^^^a = 0.05. Tc 
and T* indicate vanishing condensate density, po , and anoma- 
lous density. A, respectively. The Popov approximation leads 
to a first-order phase transition, whereas LOAF predicts a 
second-order phase transition. We have that Tc = T* in the 
Popov approximation but not in LOAF. In LOAF x' and A 
are equal for T < Tc. 



At zero temperature, LOAF predicts that the conden- 
sate fraction in the unitarity limit is 3/4, whereas in the 
Popov approximation the condensate fraction approaches 
zero asymptotically. 

Turning now to the discussion of results in the finite 
temperature regime, we note that throughout this sec- 
tion, the temperature is scaled by its noninteracting crit- 
ical value. To = {2Trh'^ /m)[p/Ci3/2)]'^^^ , where ({x) is the 
Riemann zeta function. In Fig. 3 we depict the tempera- 
ture dependence of the normal density and anomalous 
density, A, at constant p^/^a. We compare the results 
derived using the LOAF and and the Popov approxi- 
mations. For illustrative purposes, we show results for 



P 



1/3, 



1, 



,1/3, 



0.4, and p 



1/3, 



0.05. Similarly, in 



Figs. 4 and 5, we depict the temperature dependence of 
the condensate fraction, po/p, and chemical potential, /x, 
respectively, for different interaction strengths. 

We identify two special temperatures, at Tc where the 
condensate density vanishes, and at T* where the anoma- 
lous density. A, vanishes. These temperatures are the 
same in the Popov aproximation formalism, but they are 
different in the LOAF approximation. The existence of a 
temperature range, < T < T* , for which the anoma- 
lous density, A, is nonzero despite a zero condensate frac- 
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TJTn 



FIG. 4. (Color online) Condensate fraction, po/p, from the 
LOAF and Popov approximations. Similarly to Fig. 3. Be- 
cause at Tc the Popov approximation and noninteracting dis- 
persion relations are the same, the Popov approximation does 
not change Tc relative to the noninteracting case. LOAF in- 
creases Tc- 




FIG. 5. (Color online) Chemical potential, fi, from the LOAF 
and Popov approximations. Similarly to Fig. 3. 



tion, po/p, is a fundamental prediction of LOAF. In this 




FIG. 6. (Color online) LOAF predictictions for the critical 
regime, T = Tc, as a function of p^^^a: (Top panel) Rela- 
tive change in Tc with respect to the noninteracting critical 
temperature, Tq. (Bottom panel) Critical value of the normal 
and anomalous densities, x'c ~ The insets illustrate the 
p^/^a dependence of ATc/Tq = {Tc - To) /To and x'c in the 
weak-coupling regime. 



temperature range, the dispersion relation departs from 
the quadratic form predicted by the Popov approxima- 
tion for T > Tc- Above Tc the solution of the Popov- 
approximation equations becomes multivalued, indicat- 
ing that the system undergoes a first-order phase tran- 
sition at Tc- In contrast, LOAF predicts a second-order 
transition. Because at the critical temperature, Tc, in the 
Popov approximations and the noninteracting gas case, 
the dispersion relations are the same, the Popov approxi- 
mation does not change Tc relative to the noninteracting 
case. The LOAF formalism predicts a higher critical tem- 
perature than in the noninteracting case, Tc >Tq. In the 
weak coupling limit, we wave Tc — Tq, as p^^^a 0. 

As illustrated in Figs. 3, 4, and 5, the LOAF and Popov 
approximations results become qualitatively similar in 
the weak coupling limit, even though the order of the 
phase transitions remains different. However, strength- 
ening the interaction between particles in the Bosc gas 
results in enhanced differences between the LOAF and 
Popov predictions, even for temperatures, T <^ Tc- A 
larger value of p^^^a indicates stronger coupling. 

We also note for comparison purposes, that below Tc 
the large-TV approximation gives the same results as a 
non- interacting gas. Above Tc, the large- iV result above 
Tc is the same as that of the Popov approximation. Also, 
above T*, where A = in the LOAF approximation, the 
large- A^, Popov and LOAF approximation give the same 
results. 

In Fig. 6 we depict the relative change in Tc with re- 
spect to To, ATc/Tq = {Tc-To)/To, and the critical value 
of the normal and anomalous densities, x'c — ^c, pre- 
dicted by LOAF, as a function of the interaction strength 
characterized by the dimensionless parameter p^^^a. The 
insets show the weak-coupling limit of LOAF results, em- 
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phasizing the departure from the noninteracting result in 
lowest order. 

The leading-order auxiliary formalism, LOAF, pro- 
duces a more realistic set of observables away from the 
weak-coupling limit because of its non-perturbative char- 
acter. In contrast, the Popov approximation is appropri- 
ate only in the case of a weakly-interacting gas of bosons. 
The former is made explicit by studying the LOAF pre- 
diction for the relative change, /STJTq = (T^ - To)/To, 
as a function of p^^^a. The inset in the top panel in Fig. 6 
illustrates that in the weak-coupling regime, p^l^a <C 1, 
LOAF produces the same slope, 2.33, for the linear de- 
parture as that derived by Baym et aZ. [20] using the 
large-N expansion, but at next-to-leading order (i.e. they 
include density fluctuations in their calculation). The 
LOAF corrections to the critical temperature are due to 
the inclusion of self-consistent fluctuations effects in the 
mean- field x' A densities. We note that carrying 
that approach to the next order, the slope is reduced to 
1.71 [41], and is approaching the Monte Carlo estimates 
of 1.32 ± 0.02 [49-51], and 1.29 ± 0.05 [52]. It will be in- 
teresting to see how our next to leading order calculation 
compares to these results. A summary of other AT^/Tq 
theoretical predictions is found in Ref. 16. 

As the system approaches the unitarity limit, LOAF 
predicts that ATc/To 0.396 and Xc/^o AjTo 
0.873 for ^ Y. 



VIII. CONCLUSIONS 

In this paper we discussed in detail a new auxiliary- 
field formulation for the BEC problem that was first in- 
troduced in Ref. 21. At mean-field level this approach 
meets three very important criteria [16] for a satisfactory 
mean-field theory for weakly interacting bosons: (1) the 
excitation spectrum should be gapless (Goldstone the- 
orem), (2) at T = and weak coupling, it reproduces 
the known results from Bogoliubov theory, and (3) it 
has a smooth second-order phase transition. The com- 
monly used theories violate those criteria: the Hartree 
approximation violates (1), the Bogoliubov and Popov 
theories violate (3), and the T-matrix Popov theory vi- 
olates (2). Also at mean-field level, we obtain a result 
for ATc/To (Tc - %)ITq which was obtain only at 
next-to-leading order in a large-TV expansion, showing 
that including the anamolous density in our auxiliary- 
field formulation is quite important. This approach will 
be useful to study both the static and dynamic properties 
of dilute Bose gases. 

As described above, one can systematically improve 
upon the LOAF approximation discussed here by calcu- 
lating the 1-PI action order-by-order in e. The broken 
C/(l) symmetry Ward identities guarantee Goldstone's 
theorem order-by-order in e [22]. For time-dependent 
problems, however, this expansion is secular[53], and a 
further resummation is required. The latter is performed 
using the two-particle irreducible (2-PI) formalism[44. 



45]. The corresponding Schwinger-Dyson (SD) equations 
for the scalar field and the two-particle correlation func- 
tions are simplified dramatically because all vertices are 
trilinear. A practical implementation of this approach 
is the bare- vertex approximation (BVA)[54]. The BVA 
is an energy-momentum and particle-number conserving 
truncation of the SD infinite hierarchy of equations ob- 
tained by ignoring the derivatives of the self-energy, sim- 
ilarly to the Migdal's theorem[55] approach in condensed 
matter physics. The BVA proved effective in the case of 
classical and quantum A^'' field theory problems [56-58] 
and can be applied to the BEC case. In this context, 
we note that a related approximation is the 2PI-1/N ex- 
pansion which has been used in particle theory to study 
thermalization of various quantum field theories [59-61]. 
Its use for studying dilute Bose gases was discussed by 
Calzetta and Hu [38]. The 2-PI approach has been used 
also to study the quantum dynamics in the Bose-Hubbard 
model [32, 62]. 
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Appendix A: Regularization and renormalization 

Unlike the case of an operator formalism where one 
can remove vacuum energies by normal ordering, in the 
path integral method we have to subtract an infinite 
zero-point vacuum energy Vq. In addition the interac- 
tion strength Aq needs to be renormalized to obtain the 
physical scattering amplitude, as in the Bogoliubov the- 
ory for a (5-function interaction. This is accomplished 
by summing the Born series to find the physical s-wave 
scattering amplitude. We will find that regularizing by 
subtracting the leading divergences in the expression for 
the potential for the broken symmetry case is equivalent 
to dimensional regularization, which is known to preserve 
the Ward identities. It is also equivalent to renormaliz- 
ing the vacuum energy, chemical potential, and coupling 
constant. 



1. Dimensional regularization 

Our regularization scheme of subtracting the leading 
divergence is equivalent to a dimensional regularization 
procedure, which guarantees that the Ward identities of 
the unrenormalized theory are preserved. Dimensional 
regularization consists of evaluating a generalization of 
the integral in a regime where it is defined and then an- 
alytically continuing to the original ill-defined integral. 
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In the broken symmetry phase, we need to evaluate an 
integral of the form 

/[Af 2] = ^ J k^Ak ^/mk^TM^ . (Al) 

If we consider instead the integral 
1 



47r2 



dkk-" (k^ + M^y (A2) 



_ [M2]t+i r[(l-a)/2]r[(a-27-l)/2] 
- [A/2](a+i)/2 87r2r(-7) ' 

and then analytically continue this expression to a = —3 
and 7 — 1/2, we obtain the dimensionally- regularized 
value of the integral in Eq. (Al) as 



(A3) 



This is exactly what we obtained by regulating the inte- 
gral by subtracting the leading divergences, i.e. 

I^[M^]^— k^dk[^k^k^ + AP) 



because the terms we subtracted are formally zero in the 
dimensional regularization scheme. 



2. Renormalization 

In the broken symmetry phase, our regularization 
scheme of subtracting the leading divergence is also 
equivalent to renormalizing the vacuum energy, chemi- 
cal potential, and coupling constant. 

Introducing a cutoff A in the momentum integrals in 
Eq. (43), the effective potential in the broken symmetry 
case is given by 



Mo 



MoX 



2Ao cosh^ 9 Ao cosh^ 9 



(A4) 



Aosinh26' J {2n) 



and Wfe = i/cfc (e^ -I- 2%'). We first renormalize the inter- 
action strength Aq by setting 



^ d3fc 1 



Aosinh26i Asinh26i J (27r)3 4e 



(A5) 



Next we renormalize the chemical potential fiQ by setting 



Mo 



— fA6) 
Ao cosh^ 9 A cosh^ 9 ' J (27r)3 2' ^ ' 



d^k 1 



The renormalized vacuum energy is then defined by the 
equation 



Vn- 



2Ao cosh^ 9 



= Vr- 



M 



1 



2A cosh^ 9 



' d^k 
(2^ 



so that the effective potential (A4) becomes 



(A7) 



(x' + m)' 



X 
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d^fc r 1 r 
(2^1 2 I 



2 A cosh"' L 



2A sinh i 

/2. 1 



X 
2efc 



(A8) 

Inil-e-'^""']}, 



where we have taken the limit A — )■ oo since the inte- 
gral is now finite. For completeness, we note that the 
renormalized gap equation (44) for x' is 



dx' 



X^ + M 

A sinh^ 9 A cosh^ 9 



X' 



7r)3 I 



(2^)3 I 2L0k 



[2n{l3uJk) + l\ 



X 

2ujk 



(A9) 



Appendix B: Building blocks for graphs 

Mean-field perturbation theory is an expansion around 
the stationary point of the effective action and uses the 
propagators and vertices of the stationary point to con- 
struct all the graphs. The propagators that enter into the 
loop expansion are the mean-field propagators ^°f)[x]j 
where G~'^'^b[x] is given by Eq. (5), and where 
is defined by Eq. (20). The basic local ver- 
tices are the three-point vertex V*;,, which connects Xi 
with a (pa and (pb, and the two-point vertex, V^^^*", which 
changes a 4>a into a x' . The lowest-order theory also con- 
sists of the nonlocal 1-PI vertices for N-x lines, namely 



^ N 



<5^Tr[ln[g-i[x]] 
SXiiSXi2 ■ --^Xii. 



(Bl) 



These nonlocal N-x vertices are polygons made up of N 
mean- field propagators G[x]- Once we have r[(/)a,Xi] to 
some order in e, we can determine the equations for 
and X from 6T/S<pa — and ST/Sxi — s*- Subsequently, 
all higher-order 1-PI vertex functions can be obtained 
by knowing what happens when we differentiate either Q 
with respect to Xi or V with respect to both Xi and (pa- 
Because we know both and explicitly, one uses 
the identity 



SA 



-Ao 



SA- 



(5$ 



oA 



(B2) 



to obtain the rules for how to functionally differentiate Q 
and I? in a graph. Here the o symbol stands for both an 
integration and a matrix product. Using the notation of 
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Eq. (7) with Xi = Vij X'' j note that 

Functional derivatives of t/^^[x] with respect to are 
given in terms of 



5Xt 



(B4) 



where 



V\6') = cosh6i (^J , (B5a) 
V\e)^V2sinh9(^^ , (B5b) 



1 




In addition, we have 

SG[x] ^ 
^Xi 



'Q[x]oV'{0)og[x\ 



(B5c) 



(B6) 



In this notation, the inverse composite-field propagator 
2?j^^[$](x,x') defined in Eq. (20) is given by 

Pr.i [$] (a;, x') = ^ 5{x, x') + n,, x') , (B7) 

where the polarization 11*^ [<i>] is 

=0oV*^[x](6')o0 (B8) 



with 



(B9) 



= 2 [V'W°a[x]°V^(0) + V^(0)og[x]oV'(0)] . 

Another quantity we will need for obtaining the graphical 
rules is TV-^^\x\ defined by by 



oXk 

^\{^\Q)oQ\x\oV''(B)oQ\x\oV\B) 

+ v\e)og\x\ov\B)og{x\oY[())} 



(BIO) 



Similarly 



^^''\x\ 

^Xk 5x1 



(Bll) 



= i { V* o g[x] o V' o ^;[x] oV^oQoV^ 
+ V' o ^;[x] o v*^ o g o v' o g[x] o ] 



+ - { v^' o g[x] o v' o g o v'^' o g[x] o 



+ v^' o g[x] o v'^' o g[x] o v' o g[x] oV) . 



We also define the leading-order 3-x 1-PI vertex function, 
QT[x]. as 



ciXfe oxk 

= -<t>oTv^'^[xmo<f> 



(B12) 



+ - Tr[g[x] o v\e) o g[x] o v^(0) o g[x] o v^{e) 
+ g\x\o v\6) og[x]o v^{6) og[x]o v\e) ] . 

The 4-x vertex is then given by 

= cj, o TV^'^^xm ° - I Tr[g o V*((?) (B13) 

og oV'{e)og o v'' {e)ogo v' {e) ] + perms. 

With the above definitions we can construct the rules 
for inserting a or x vertex into a graph: Inserting a x 
line into g[x\i we obtain: 

= ^g\{x,.z,) v-m g\iz,,x,) 

0Xt(zij 

(B14) 

= -gov\e)og. 

Inserting a x line into 2?[<I'], we obtain 

SV[x^Y^^{z,,Z2) 



5Xk{z3) 



(B15) 



[dz4] [dz5]2?^"(^l,^4) ^%^^^^^2?"-'>5,^2) 



SXkizs) 



We also need to insert a line into V. The 2-x l-i/i vertex 
is given by 

^3„. , ^ SV-^^^{z, z,) ^ V^%iz,,z,) riz2) 

+ 0=(Z1)V^^%(Z1,Z2)5(Z2,X1), 

and for the 2-Xi '2-4>a vertex we find 

r^^:,at ^ S'V-^^HZ„Z,) ^g^g^ 
S(f>a{xi) 5(i)b{x2) 
^ 5{ZI,XI)V'^''\Z^,Z2)6{Z2,X2) 

+ (5(zi,x2)V^J"^(zi,z2)<5(z2,a;i). 



Thus we obtain 

5V''i[mzuZ2) 



(BIT) 



5(t)a{xi) 

= -J [dz3] [dz4] T^'"'{Z1,Z3) r^„'^(z3, ^4, Xi) I?"^'(Z4, ^2) 
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1. Inverse propagators to order e 

Using the above rules, we derive the one and two-point 
vertex function to order e. For the one-point function in 
the presence of sources we have the following two equa- 
tions: For (h't we have 



(B18) 



-Tr[2?r 



3 ..a ■ 



whereas for Xi fi^id. 

s = —f. = X 

r . 1 



o V o , 



(B19) 



, , „.Tr[goV'] + -Tr[I?og3-] 



o- 



r = — ^ + 



rrt= Y [(<|)°V) + (Vo,^)i 



$3 



Q3>- 



Q4 



+ (9(8^) 













r4 







+ O(e-) 



FIG. 7. Feynman diagrams for the second derivatives of F. 
Solid and wavy lines correspond to the propagators of </!> and x- 
Dashed lines denote (h. 



In turn, for the inverse propagator matrix we have 

52r[$] 



(B20) 



and 



-lab 



(B21) 



- ^Tr[L'oF^- oDoF;^-] + ^Tr[DoFf-] 



e 

2i 



The term that mixes (p and x is 



(B22) 



^Tr[P o o I? o F-- ] - ^Tr[I? o [i?- o 0]« ] 



The propagators for the theory are obtained by inverting 
the 5x5 inverse propagator matrix. Expanding the prop- 
agators in a power series in e and keeping terms to order e 
gives the graphs for the propagators that one would have 
obtained by working directly with InZ to order e. The 
Feynman diagrams for the second derivatives of F[<i>] are 
shown in Fig. 7. 



2. n'^[$](x,s') 

To conclude this section, we complete the calculation 
of n'-' [(f>](x,a;') introduced first in Eq. (B7) and explic- 
itly evaluated in Eq. (B8) above. In the imaginary-time 



formalism, we first introduce 

M'^[x]{x,x') 
1 



(B23) 



TY[g[x]{x',x) o v'{e) o g[x]{x,x') o v^{9) ] 

6 



(27r)6 



ni ,712 — — oo 



1 



X - Tr [ ^ [X] (k2 , na) o V* (0) o g [x] (ki , m ) o (9) [ 

X g«[(kl-k2)-(r-r')-(w„j-w„2)(r-T')] 
L I a k ^ gi[k-(r-r')-w„(T-r')] 



/? J (27r)3 ^ 

^ -' n— — oo 

1 /"/"d^fcid^fca 



/3yy (27r)^ 



^ (2^)3^(k-ki-f k2)(5„^„,. 



1 



X - Tr[g[x](k2, n2) o ^(9) o g[x](ki,ni) o V^(0) ] . 

Expanding Ai^^ [xli^i in ^ Fourier series, 
M'^[x]{x,x') (B24) 



1 f d^k 



J2 M'^IXU^, n) e«[k-(r-r')-<^„(T-r' 



)] 



where from (B23), TVI*^ [x](k, n) is given by the convolu- 
tion integral. 



A^-[x](k,n) = i, / ^ Y: 



(B25) 



2(3 J (27r)3 ^ 

n' — — OO 

X Tr[^[x](k - k',n - n') o V^^) ° ^[x](k',n') o V^(0) ] 
From Eqs. (B8), (B9), and (B23), we then have 
ff^[$](k,n) = 0oV^J'[x](k,n)o0-7W'J'[x](k,n), (B26) 
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with 

V^^ [x](k,n) = ^ [V^(0)og[x](k,n)o V^"(0) (B27) 

+ v\e)og[x]i-k,-n)oV'ie)] . 
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